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Abstract 

The theory and experiments concerned with the electron-ion thermal relaxation and melting of overheated crystal lattice constitute 
the subject of this paper. The physical model includes two-temperature equation of state, many-body interatomic potential, the 
electron-ion energy exchange, electron thermal conductivity, and optical properties of solid, liquid, and two phase solid-liquid 
mixture. Two-temperature hydrodynamics and molecular dynamics codes are used. An experimental setup with pump-probe 
technique is used to follow evolution of an irradiated target with a short time step 100 fs between the probe femtosecond laser 
pulses. Accuracy of measurements of reflection coefficient and phase of reflected probe light are 1% and ~ 1 nm, respectively. 
■It is found that, firstly, the electron-electron collisions make a minor contribution to a light absorbtion in solid Al at moderate 
intensities; secondly, the phase shift of a reflected probe results from heating of ion subsystem and kinetics of melting of Al crystal 
during < t < 4 ps, where t is time delay between the pump and probe pulses measured from the maximum of the pump; thirdly 
the optical response of Au to a pump shows a marked contrast to that of Al on account of excitation of of-electrons. 

Key words: femtosecond laser ablation, pump-probe, optics of hot Al and Au 
PACS: 52.38.Mf, 52.25.Os, 02.70.Ns 
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s^i. Supersonic heating and melting 

Figures II 121 show diagrams of processes in pump femtosec- 
ond laser pulse (fsLP) action on metal. The three time slices 
I "ei", mim2, and C1C2 in Fig. Q] correspond to the following 
non-equilibrium processes: (e-i) the electron-ion thermal relax- 
ation, (m) the melting of an overheated crystal lattice, and (c) 
the cavitation decay of a metastable state. Duration of fsLP 
Tl ~ 40 - 100 fs is shorter than characteristic times of these 
three processes. They have very various time scales from subpi- 
coseconds to nanoseconds. The electron overheating (Te » T, 
starts from eij when a fsLP arrives J3, 0, 0. 0. H S @ 
and disappears at ei2 when temperatures T e , T, equilibrate 
(teg — t e a - 3 - 6 ps for Al at our intensities). The time is 
• *h reckoned from the maximum of pump fsLP in Fig.Q] Since ar- 
/\ riving of the pump to a target the conductivity electrons become 
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much hotter than the ions. 

Two-temperature (2T) matter with hot electrons transits to 
a peculiar state with thermodynamic and optical characteris- 
tics different from one-temperature (IT) case. In 2T there are 
appearance of excesses of electron energy and pressure above 
equilibrium IT ones. Also there are changes in elastic moduli 
and band structure. In semiconductor lattice the binding forces 
become weaker with increase of T e , while in metals situation is 
opposite. Large changes in optics of Au at high T e result from 
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Figure 1 : The pump fsLP tz, and a chain of kinetic or transient processes 
(l)"ei", (2) "mim2", and (3) "C1C2" initiated by absorption of the pump. 



excitation of of-electrons. On account of the ion heat capac- 
ity C, (thermal "inertia" of a lattice) the beginning of melting 
t m \ ~ CjTm/aTe is delayed relative to the instant "eii", where 
T m and a are the melting temperature and e-i energy exchange 
rate. 

It has been known that: 

(i) for metals and semiconductors the fluences F near the 
ablation threshold Fm are significantly higher than the melt- 
ing threshold F m lllOj], In this sense the Lennard-Jones (LJ) 
case is an exception 1 1 Oil - LJ E m /E co h vacancy migration to 
cohesion energy ratio and T-}/T c triple to critical temperature 
ratio are large in comparison with metals and semiconductors. 
While Eabi/Ecoh is approximately the same for all three groups. 
Therefore it is not suprising that in the LJ case F a t,i ~ F me [ t . 
In metals for F > F a u the molten layer is approximately as 
thick as heated dj. When stretching stress overcome mate- 
rial strength a cavitation (fragmentation of liquid) begins inside 
molten metal at the instant c\ shown in Fig. [TJ In LJ near 
threshold F ~ F„bi spallation (fragmentation of solid) starts in 
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Figure 2: Electrons are heated up by a pump in a skin layer of metal (the rect- 
angular "skin"). The "EHC" supersonically (comp. the heat dj and sonic 
propagations c s t e „ at t = t eq ) carries heat from the "skin" layer into the bulk 
of the target along the horizontal arrow (it shows the heat flux) forming the 
heated layer dj . The energy of hot electrons is transferred to ions in fast pro- 
cess of "e-i heating" (the vertical arrows) [3]. As a result the crystal lattice 
is overheated to AT above the solidus temperature. The solidus and liquidus 
temperatures at p° are 1.2 kK, 1.48 kK (Al), 1.9 kK, 2.05 kK (Au). They are 
significantly higher than the melting or the triple point temperatures 933 K (Al), 
1337 K (Au) corresponding to melting at low pressure. 



deformed crystal f^fo. Above threshold F > F a bi the LJ spalla- 
tion transforms to cavitation as molten LJ layer becomes thicker 
and the fragmentation zone transits from solid to molten LJ. 

(ii) the electronic heat conduction wave "EHC" in Fig. [2] is 
supersonic within the 2T slice "ei" shown in Fig. [T] 0,13]. 

These facts (i,ii) result in the isochoric heating and stress con- 
finement [2, 6]. Estimates of the EHC speed are: xehc ~ yjxi, 
X = /v/3 is a thermal diffusivity, I = v/v ~ 1 nm is a mean 
free path, v = 1/t is a collision frequency, v is the Fermi 
velocity. Therefore the Mach number of EHC wave is high 
xehc/ c s ~ 100 Vr/f up to a few picosecond. Within the time 
period ~ t eq the "EHC" creates a heated layer dj « 100 nm 
thick in Al and * 250 nm thick in Au H 1 Q % 0, % . 

The dj is much thicker than acoustic penetration depth c s t 



at t < teq as illustrated in Fig. |2] As a result of (i) and (ii) 



for a force balance, and two thermal equations 
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describing instant local electron and ion heat balances, where x 
is an Eulerian coordinate defined in Fig. 12 x(x°, t) is the trajec- 
tory of material particle with Lagrangian coordinate x° = x(t = 
—5tl) which is equal to x before the action of a pump fsLP. 
In 2T-HD simulation the pump intensity is I{t) = 0, t < -3tl, 
and /(f) = [F/(t l V^)]exp(-f 2 /Y£), t > -3r L ; p° = p(x,t = 
-3tl) is the initial material density: 2.71 g/cm 3 (Al) and 19.3 
(Au), p,p are pressure and density, k is a coefficient of a ther- 
mal heat conduction. Q(x,t) = [I(t)/5] exp(-x/6), where 5 
is the thickness of the skin-layer shown in Fig. [2] as "skin". 
Unknowns p, T e , Tj, u are functions of the Lagrangian vari- 
ables x°, t. When knowing the trajectories x(x°, t) , the vari- 
ables p, T e , T,, u can be presented as the functions of the Eu- 
lerian coordinates x, t; E e ,Ej are the internal electron and ion 
energies, Q is the absorbed power, a{T e -Ti) is the e-i energy 
exchange term. 

From the profiles p, /, T e , T,, as functions of x , obtained 
from the Eqs.(l-3), one can find the corresponding profiles of 
the dielectric permittivity s(p, /, T e , 7 1 ,) - s r + /e, and the com- 
plex index of refraction N - n + ik, s = N 2 , where / is 
a volume fraction of liquid phase discussed below for case of 
solid-melt mixture. Next, for the x -profile s(x,tfi X ) at fixed 
time tfi x (velocity of light is taken to be infinite), Helmholtz 
equation 

(4) 



d 2 F/dx 2 + ek 2 F = k = co/c 



is solved for amplitude F of the probe fsLP perpendicular to 
a target. Equation (0]) describes a reflection of the probe light 
from target. Its solution gives the amplitude and phase of the re- 
flected wave and correspondingly the time evolution R(t), \fj{t) 
of the reflectivity R and phase They are compared with 
the experimental dependencies in Figs l3l4l By contrast, the 
pump absorption is taken from the experiment J5t]. Equation 
(0]) is evaluated by the 2x2 transfer matrices method I19II . In 
our experiments the chromium-forsterite laser with the pump 



there is a volume non-equilibrium melting in the slice "miir^" 
J2l with formation of overheated solid grains surrounded by 
melt. At the slice "m 2 - fr" t eq < t < "fr" in Fig. [T]the A pump = 1240 nm and probe A pro b = 620 nm (first and sec- 
ond harmonics) is used A fsLP duration is 100 fs. 
Values R, if/ have been measured by microinterferometric tech- 
nique described in H HE H3]. 

3. 2T thermodynamics, collisions, and thermal transport 

The values p,, £,(p, T,) and p e , E e (p, T e ) in (1-3) are taken 
as in ^ from the wide range equation of state j2lll and from 
Fermi model for the conduction electrons, p = p, + p e is the 
total pressure. The coupling factor a and the heat capacity C e 
are taken from |[§] . According to the Drude formula the electron 
heat conductivity k in (2) is 



heat wave velocity becomes much lower than c s - and then the 
well-defined melting/recrystallization front is formed |2ll6i 1 1 ill . 
At t — "fr" ~ 1 ns the melt layer covering a residue lll2ll of the 
target is completely solidified. 



2. 2T hydrodynamics and light reflection 

In Lagrangian variables the 2T hydrodynamics (2T-HD) 
equations are 

dx{x°, t)/dt = u, pdx(x°, t) = p°dx" 

for kinematics and mass conservation, 

p°du/dt = -dpldx", (1) 



*r=(l/3)v*C e /v, v = {v d 2 g + v-f) 



(5) 
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Figure 3: Role of the e-e collision frequency v ee in optics. The solid curve 
obtained for v ee = 0. The experimental (squares) and theoretical reflectivity 
R(t) normalized to initial Rq ofAl; Z = 3, m e ff/m e = 1.5, A/ ? t, from flUl . 



v = + 3kBT e /m e , \f = y/2Ef/m e , Ef = ksTf is Fermi 
energy. At T e < Tf when electron degeneracy is significant the 
collision frequency is v « Vd eg ("deg" stands for degenerate), 



Vdeg = (v ei + v ee ){plp°) 



■1.3 



(6) 



The factor (jo/p ) -1 3 in (O approximates the quantum- 
mechanical molecular dynamic (MD) data 12211 showing the 
drop of the electrical conductivity of Al with the density de- 
crease in the temperature range under consideration. The 
electron-electron collision frequency is 



v ee = b(Ef/h)(T e /Tf) 



(7) 



The Coulomb collisions v p i in (0 is taken as v p i = 
(E F /n)(T e /T F )- 3/2 (p°/p) 2/3 . They dominate at very high T e > 
Tf and limit v (0 at T ~ by an atomic frequency 
v„, ~ 10 16 s~' (saturation of the v increases with increasing 
T e ). 

The electron-ion collisions v e , in va eg 
separately for solid and liquid Al 



5T6J is calculated 



Jiq 



4.2 • io 14 (7yr 3 ) (8) 

1.1 ■ 10 14 r,/(130 + 0.03677,- - 66700/r,) [ s -\ (9) 



where T$ = T,„(p — 0) = 933 K, T, in K. The coefficients 
in (0 approximate the quantum-mechanical MD heat conduc- 
tion of the molten IT Al 112 311 . The thermal conductivity is 
calculated with Kubo-Greenwood formula up to T = 10 kK 
1 2311 . Values for k from (151618b and for electric conductivity 
cr = n e e 2 /(m e v e j) with ([8]l describe well IT reference data for 
solid Al between the Debye and triple point T3 temperatures. 
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Figure 4: The drop of the phase shift Ai/< results from the solid-liquid transition 
within first few picoseconds, Ai/r = if/(t) - if/Q, if/Q is an initial phase (a phase 
of reflected wave from a cold target). Diamonds indicate experimental data. 



4. Optical properties: collisions and interband excitations 

The value e necessary for (O is a sum of the Drude and the 
interband terms [l7l 0, 0] 
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where 



co , = Ann e e /m e f f , 



in Al m e ff = (1.2— 1.7)m e . In an Al crystal results mainly 
from the transitions between parallel zones 111 811 . This term 
dominates at room temperatures. Its contribution increases the 
absorption more than by order of magnitude in comparison with 
pure Drude absorption. In a molten Al A^, disappears 12511 . 
Electron density of states (DOS) in Al is stable against melting 
II26J] and against T e increase as was checked up to T e = 70 kK 
in 12711 . Therefore the ion charge Z = «<>/«; and the electron ef- 
fective mass m e ff defining a> p i weakly depend on melting and 
overheating of electrons. A phenomenological dependence of 
Afei (TTOTb on a total frequency v i5\6\ has been proposed 
At large v > co the term A^ becomes small as compared to 
the Drude term as it is in liquid. 

5. Optics of mixtures 

MD simulations show that in Al within the early (~ 0.1 - 
few ps) stage the solid-liquid mixture fragmentation space scale 
~ 1 nm is small in comparison with 8. Therefore the s" ux 
can be defined by a volume fraction / of liquid in mixture: 
s"" x (f, s so1 , e llq ). For the weak mixtures (/ « or / » 1) we 
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An approximate interpolation for intermediate / 
linear in / solutions (fTTT i is 
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1 of these 



/ 4 ). 
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Expressions (11 U12t have been used in calculations shown in 
Figs.[3]and|4] At an early stage it is necessary to consider optics 
of a solid-liquid mixture given by these expressions because at 
this stage the thickness of a solid-liquid mixture layer is compa- 
rable with the thickness 8 of a skin-layer - a penetration depth 
of a probe photon. 



ImF 



Values s sol (T = 300 K) = -53.5 + 24. li, N so1 = 1.6 + 7.5z, 



sol 



E n i(T = 1200 K) = -40 + 15.3/, N liq = 1.2 + 6Ai 125 [differ 
moderately. In this case there is another approximation 12811 



rnix\ 1/3 



(s m,x ) 



,1/3 



nonlinear in / but linear in small difference s so1 - s bq . Com- 
parison of this approximation with (11 1112b shows that results 
differ small: the maximum deviation achieved at / « 0.75 is 
« 1% for Re (e) and less than 3% for Im (e). As was said at 
early stage optics of solid-liquid mixture is significant. Let us 
mention that as it will be shown below at this stage there is non- 
equilibrium melting of an overheated crystal. At intermediate 
and late stage C1C2 shown in Fig. Q] existence of vapor-liquid 
mixture can influence optical reflection if thickness of cavita- 
tion layer JUI^O is less or comparable with 6. The cavitation 
layer covers the undersurface vapor-liquid layer against probe 
photons. 

6. Role of v ee in optics of aluminum 

In an early stage the T e is high - then the v ee (0 dominates 
in the total frequency v |5) as maximum T, for our fluences 
are smaller than 10 kK and the v e , collisions are less frequent. 
There are different coefficients b K and b s for the v ee (0 in 
expressions for k (O and for optics - in the Drude and in the 
interband terms (TTOt . In crystals the b K includes normal and 
Umklapp processes while the b s in the solid Al may differ from 
zero as a result of the Umklapp effect (situation is different for 
gold, see below). In liquid Al the coefficient b s = - the 
e-e collisions do not contribute into optical absorption as the 
Umklapp is impossible. In Al a Fermi sphere is larger than in 
Au while the Brillouin zones are approximately equal (the lat- 
tice constants are « 4 nm for Al and Au both). Therefore the 
Umklapp effect is more significant in Al. E.g., it results in or- 
der of magnitude increase of the electron-ion energy exchange 
rate a 112911 . Value of b in (0 is a subject of discussions. Be- 
low the Debye temperature the v ei oc T 5 tends to zero strongly 
and becomes less than the v ee oc T 2 . Then specific electrical 
resistance r — l/cr for very pure crystals is r — AT 2 + CT 5 . 
The measurements |30i| give b e = 15. Here we suppose that 
v ee does not depend on frequency of electromagnetic field and 
the coefficient Z? in (|7j for resistance r is equal to b s . At the 
same time the theory lf30h gives b B = 0.6. For high T e the 
calculations |3l b E — b" — 1 is accepted. 

Estimate of the upper limit for b K follows from the check of 
the Wiedemann-Franz law k/ct = LT done for melt Al in l23ll . 
where L is Lorentz number. In Drude approximation we have 
k/ct = LT/(l + v^ e /v e i) because in melt umklapp is absent and 
v ee (7) does not contribute to electric conductivity cr = 1/r. 
Relative deviation e of the ratio tc/ (oT) from L in 12311 is less 
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Figure 5: Deviation of a phase of reflected light resulting from melting, F is 
a complex amplitude of a reflected wave. = n - 12.5 [nm]47T/Ap mb and 

- 14.7 [nm\AnlAprob- 



than 10%. This means that electrons in Al remain degenerate up 
to temperature T = lOkK achieved in 1E3I1 . The limit e < 0.1 
impose restriction on value b K < 1.5-2. It is obtained from 
v£ e /v' ! ? < e where expressions (7) and (9) was used. On the 
other hand, the value b" is important at early stage because it 
influences heat propagation into bulk when the propagation is 
supersonic. To achieve thickness dj « HOnm at Al ablation 
threshold the values b" should be near this restriction b K « 
1.5-2. The thickness dj « HOnm is necessary to reproduce 
experimentally defined crater depth 45 - 50 nm. For smaller 
b K the simulated crater is deeper and fluence threshold is higher 
than the experimental ones. 

Performed optical measurements shown in Figs. [3]|4]together 
with simulations give additional information about b e . The in- 
teresting sharp narrow "well" at the theoretical R(t) in Fig. [3] 
corresponds to b s = 1.7. Its minimum is achieved when T e 
and hence v ee [7] have the largest amplitudes (at the end of the 
pump). This well might be very useful for diagnostics. The left 
wing of the well follows the history of electron heating while 
the right wing reflects the kinetics of melting because the grad- 
ual phase transformation from crystal state to melt in a skin- 
layer progressively suppresses optical contribution of the e-e 
collisions (0. Unfortunately the well is not observed. Density 
of experimental points at the time axis in Figs. [3] IU is large 
enough to exclude missing of the well between the two succes- 
sive points. The accuracy of experimental measurements of a 
relative reflection (« 1%) is sufficient to catch the well. Anal- 
ysis of our simulation runs in a fluence range 1 < FJ Fm < 4 
shows that it is necessary to have the b s below 0.2-0.5 to 
meet the measurements; for Al the calculated and measured 
thermomechanical ablation threshold is F a i,i\i„ c - 0.75 J/cm 2 , 
FabLbs = 65 mJ/cm 2 J3]. 

7. Melting and decrease of phase angle 

Simulations show that the phase evolution i(/(t) presented in 
Fig. [4] contains important information concerning the kinetics 
of melting. The base for this is the difference commented in 
J23I1 . Fig. [5]between the defined in Section 4 e™' and e'' 9 . As a re- 
sult of attenuation of the band-band transition during melting 
the value 11 = Re(N) becomes smaller. This is why the i//l q in 
Al is 2.2 nm differing from solid in the direction of the phase 
rotation to the crater. The value 2.2 nm corresponds to the case 
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Figure 6: Rise of 7*,- in 2T-HD and MD (the fluctuating profiles) simulations, 
Al, Fabg = 65 mJ/cm 2 . 
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Figure 7: Evolution of the phase composition as a result of the increase of Tj. 
The e-i heating gradually rises T, as it is shown in previous Figure. Transver- 
sally averaged .s -profiles obtained from MD simulations are presented. The 
arrows "liquidus" and "solidus" at t = 4 ps are taken from previous Figure. 
The values of the MD symmetry index s(x, t) = 2.5 - 2.7 taken at the 2T-HD 
solidus fronts are approximately the same for different instants. 



of a Fresnel reflection from homogeneous semispace. Remark- 
ably that this small difference is measurable by the pump-probe 
interferometry. The sign of rotation directions in Fig. [5] de- 
fines the sign of the phase difference AiA_in Fig. |U Detection 
of ultrafast melting of semiconductors 13 If is possible due to 
the same liquid-solid phase difference At//. But in this case the 
melting transforms semiconductor into a metallic state - there- 
fore At// is significantly larger (e.g., At// = 12.4 nm for GaAs, 
A = 620 nm) - and this transformation can be detected easily. 

The phase At//(t) obtained from 2T-HD equations is com- 
pared with experimental data in Fig. [4] Expansion movement 
of reflecting boundary should increase At// but at the early time 
it decreases as a result of gradual melting of skin-layer. Agree- 
ment between data and theory indicates that theory given below 
properly describes the melting. 

Figure|6]illustrates the heating of ions by hot electrons in 2T- 
HD model. In MD simulation atoms are heated by the space- 
time distributed thermostat power source with the temperature 
distribution taken from 2T-HD. This MD approach is similar 
to one developed in J2, I32I1 . see also where phenomeno- 
logical terms describing non-equilibrium melting are added to 
2T-HD equations. Particular mechanism of heating (heat flow 
from electrons or thermostat) has no action upon the kinetics 
of melting if we suppose that elastic moduli do not depend on 
T e as in the case of Al |27|. In Fig. the 2T-HD and MD T t 
profiles are approximately the same. Some difference results 
from equilibrium and non-equilibrium description of melting. 
The arrows in Fig.[6]mark the slice of melting from equilibrium 
2T-HD. The kinks at the ends of this slice result from hidden 
fusion energy. In MD the degradation of crystal symmetry dur- 
ing the fast heating and melting is distributed in wider range 
beyond the liquidus/solidus positions. Crystal beyond solidus 
is in overheated state ifll^. 



Phase transformation and propagation of melting into bulk 
is shown in Fig. [7] The profiles of the symmetry index s are 
presented. The index s is defined as a number of crystal axis 
passing through an atom and averaged over atoms within a x- 
slab, and 5 = 6 in a cold fee lattice. One can see how quickly 
the rather thick (thicker than S) layer of mixture is formed. 
Later the layer of pure melt (the plateau at the instant t — 4 ps 
s— profile in Fig. [7|i appears. Much later the narrow melting 
front with small overheating separating a melt from a crystal is 
formed. Maximum thickness of the molten layer for the fluence 
F a bs - 65mJ/cm 2 is « 100 nm. 

8. Gold optical response and rf-electrons 

Absorption of pump fsLP causes sharp changes in optical 
properties at very early time - during the pulse T/,. They are 
shown in Fig. [8] Bulk gold targets are used. This response 
is caused by fast heating of electrons. If we compare Figs. 
13 141 (Al) and Fig. [8] (Au) having similar relative temperatures 
T e /Tp at 2T stage and similar final T,- = T e temperatures af- 
ter e-i relaxation we will see obvious large differences. They 
are related to the differences in the band structures of Al and 



Au 12711 . Estimates of heating history from the maximum T e 
to e-i thermalization give for the main seven thermal param- 
eters F 'ine j F \nc\ablt F[ nc J/cm 2 , F ahs mJ/cm 2 , E e \max MJ/kg, 
T e \ max eV, Z = N e6sp , Ti\ max kK the values: (0.5, 0.7, 50, 1, 
1.5, 1.2, 0.8 ), (1, 1.3, 140, 3, 2.5, 1.5, 2.5), (2, 2.9, 500, 10, 5, 
2.4, 8) for the three cases shown in Figs l8l9l where Fj„ c \ a i,i = 
(1.3 - 1.4) J/cm 2 is a thermomechanical ablation threshold, 
E e \max> T e \ max are electron thermal energy and temperature at 
the end of a pump, Z is number of electrons exited from the 5d 
to the 6s,6p band at the maximum T e \ max , Ti\ max is maximum 
ion temperature achieved after e-i relaxation. An expression 
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Figure 8: Sharp change of R and ip during pump action. 

E e (T e ) = 45.7499r ( Ml9.756r 21 + 105.419r 2 - 2 -30.955ir 2 - 3 
for electron thermal energy at fixed density p - 19.3 g/cc ap- 
proximates data obtained from the ABINIT [33] simulations up 
to T e = 10 eV, here E e is in MJ/kg, and T e in K. Our curve 
is above the standard parabola E e = because we in- 

clude exciting of the cZ-electrons. But it is below the curve |Ql 
where a red-shift of the c/-band with T e was not considered. 
Experimental and theoretical values for the depth of a crater at 
a threshold F;„ c |aM is 1 10 nm fl7[] . The first case with the small- 
est Fjnc is near a melting threshold for bulk Au. Our three cases 
cover a range of energy densities obtained in 1 1_5 , 34 , 35 , 3(| 37 ] 
for ultrathin (25-30 nm) freestanding Au films. Here we con- 
sider bulk targets and use different technique of measurements. 

The values of Z presented above have been defined using the 
ABINIT 13311 code, the normalization condition for the number 
of electrons, and the expression for the amount of electrons in 
the 6sp band. The ABINIT has been used as in l27ll for calcu- 
lation of the T e dependent DOS at different T e supposing that 
density is equal to p° (isochoric heating). For given T e the 
DOS of the 5d band obtained from ABINIT has been approxi- 
mated by a rectangular and the DOS of the 6s,6p band has been 
described as the function g(e) oc -\fi in order to calculate Z. 
At a given T e a root p of the normalization condition 



11 = 



Vim 3 / 2 (k B T e )^ 

7T 2 n a tnm 



f 

Jo 



yfxdx 



e x-n/k B T e + 1 
,(ji-E l )lk B T e 



+ L, 



(13) 



where L = (10/C^i -E 2 )) ln((l +e ^-^» (I * B '«))> 

defines a chemical potential p(T e ), here E\,E% are edges of 
the 5d band relative to the bottom point of the 6s,p bands. The 
first and the second terms correspond to the numbers of elec- 
trons in the 6s, p and 5d bands, resp. At room temperature these 
numbers are Z = 1 and 10. It is known that for Au differ- 
ence between this approach and calculation of p(Te) with ex- 
act DOS is small 113 811 . The exact function p(T e ) is obtained 
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Figure 9: Fast growth of £ as result of electron heating. The £ = £,- + isj, £, < 
0, e, > was calculated from data shown in Fig. [5] and the Fresnel formulae 
valid while the £ -profiles may be approximated by a step function. The solid 
curves presents data from 1391 at E e = 2.9 MJ/kg. They corresponds to our 
intermediate case Fj„ c = Fabl< Two bars give £ for liquid Au at 1337 k |2^1 . 



together with the DOS in the ABINIT simulation. The excita- 
tion degree Z is given by an expression Z(T e ) = 11-L, where 
L stands for the second term in ( TTSl but now with known p(T e ). 
To estimate possible influence on the value of Z of position of 
the bottom point of the 6sp band this position was varied in the 
range +2 eV relative to its position at T e = 0. 

The values of Z and collision frequency v are necessary for 
the Drude estimates of s. Growth of them is responsible for 
the rise of in Fig. [9] At room temperature v/u) pro b is 

small: 1.2% from electrical and thermal conductivities or 3.3% 
from optical data J39tl . ha> pro b — 2 eV. Heating of an electron 
subsystem in our conditions rises v to large values: v/cj pro t, ~ 
1 . There are three candidates responsible for the growth of v : 
(i) an enhancement of v e , in solid or liquid Au with T e » T, 
in comparison with T e = T, case, (ii) v ee and Umklapps in 
crystal Au, and (iii) the e-e collisions between electrons from 
the d and the sp bands in solid or liquid Au. Let's consider 
these cases. 

(i) The frequency v e/ and the e-i exchange rate a are cou- 
pled, Eej = a(T e - T t ) » aT e ~ n e AE ei v ei (T e ), T e » 
Tj 14011 . where AE e i ~ (m e /Mj)Ef is an energy trans- 
fer from a light electron to a heavy ion in one e-i colli- 
sion, v e j(T) ~ a 2 (T/9) yfmjWi, a is an interatomic dis- 
tance in condensed matter, 6 is Debye temperature. There- 
fore v e i(Tj) ~ aTj/(n e AE e j). If a — const as in Al ifsL l29ll . 
then v e i = v e j{Tj) oc Tj as usual. But in gold there is sig- 
nificant dependence on T e , a — a(T e ). An approximation of 
calculation of a oc 

/ (g(e)) 2 (-df/de)de with DOS g(e) from 
ABINIT Q (including the shift of the d-band with T e ) is 
a(Te) = (0.23 + 4.3r e 3 6 /(l + + 0.9r f 41 ))10 17 W/m 3 /K. 
It is valid up to T e = 10 eV, here /(e) is Fermi distribu- 
tion The a(T e ) begins to grow at T e — 3 kK from 
a(T e < 3 kK) = 0.23 ■ 10 17 W/m 3 /K and saturates at the 
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value ten times larger » 2 ■ 10 17 W/m 3 /K 14 ill . This means 
that v e j{Ti,T e ) ~ a(T e )Tj/(n e AE e j) is enhanced if T f » T/, 
r c > 3 kK. There is v s °'(T) = 1.2 • 10 U T s"\ v'^(T) = 
3.3 • 10 14 + 1.5 • 10"r in IT gold, T in K. For T t = 2kK we 
have v"? = 6.3 • 10 14 s -1 . \f T e > 2 eV as in two cases shown 
in Figs l8l9l then the enhancement due to increase of T e gives 
large frequency v e i(Ti = 2 kK, T e — 2 eV) ~ 6 • 10 15 ~ 2u) pr0 b. 

(ii) In Al v'"" kl is small, Sec. 5. As was said, the Umk- 
lapp processes are weaker in Au as the result of smaller 
Fermi/Brillouin ratio. Therefore it seems plausible that this can- 
didate is less important than two other candidates. 

(iii) Electrons from d versus sp bands have different angular 
momentum and effective mass. Therefore photons are absorbed 
in collisions between them. Corresponding frequency is 
given by (0) with possible dependence b(T e ). The may be 

~ Wprob at T e ~ Tf. 

From imaginary and real parts of expression dTOb we ob- 
tain v = v/a) prob = (e,- - A,)/[l - (s r - A,.)] and Z/m eff = 
(1 + v 2 )(e i - A / )/20 j 6/v. The mass m^ f /m e = 0.95 - 1.15 at 
room temperatures [39] remains approximately the same with 
increase of T e . This follows from our ABINIT simulations. 
According to l36ll in 2T Au Ed > ti(jJ pro b — 2 eV, where Ed 
is an absorption edge of the li-band. If we neglect the band- 
band term A h i, = A,. + ;A ; in ( [Tol l. A r - 0, A,- = 0, then 
(Z; v) = (0.9; 0.5), (1.7; 0.9), (3.3; 1.3) at the maximum \s r \,Si 
in the three cases shown in Fig. [9] 

Our data shown in Fig. [9] agree with data from l36ll in e, 
but give larger e, and have a maximum at the time dependence 
\e r (t)\ while the dependence e r {t) from J36ll saturates. Perhaps 
the last difference results from the conductive cooling absent in 
ultrathin films. Let's mention that data l36Tl contain the useful 
dependence s(ai) but may be less accurate at a particular fre- 
quency. (Z; v) = (1.2; 0.4) for the maximum of the dependence 
e(t) from 1 36] shown in Fig.|9]as the solid curves. As a result 
of smaller e, these values are below than our values (1.7; 0.9) 
for the case with approximately the same energy E e . Never- 
theless there are appreciable excitation Z and frequent colli- 
sions. Therefore we can conclude that measurements confirm 
the theoretical findings presented above that a pump irradia- 
tion creates an excited population (Z > 1) rising <d p i ( TTOb and 
transfers gold into the state with strongly collisional widened 
energy levels. 
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